Crossover physics in the non-equilibrium dynamics of quenched quantum impurity 

systems 
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A general framework is proposed to tackle analytically local quantum quenches in integrable 
impurity systems, combining a mapping onto a boundary problem with the form factor approach to 
boundary-condition-changing operators introduced in Phys. Rev. Lett. 80,4370 (1998). We discuss 
how to compute exactly two central quantities of interest: the Loschmidt echo and the distribution 
of the work done during the quantum quench. Our results display an interesting crossover physics 
characterized by the energy scale TJ, of the impurity corresponding to the Kondo temperature. 
We discuss in detail the non-interacting case as a paradigm and benchmark for more complicated 
integrable impurity models, and check our results using numerical methods. 
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Introduction. Quenches are a clear-cut way to study 
the non-equilibrium physics of quantum systems and re- 
veal their intrinsic time scales. Here, energy is injected 
into an otherwise closed system at time t = via a sudden 
change of a control parameter, and it is subsequently left 
to evolve unitarily. From a theoretical perspective this is 
convenient because energy exchange between the system 
and its environment is not considered during its evolu- 
tion. Furthermore, such quench scenarios can be realized 
by experiments probing the transient dynamics of physi- 
cal systems, such as atoms in a tunable trap potential or 
quantum dots connected to tunable control gates. 

In many cases, the effect of interactions is to give rise 
to crossovers between two different kinds of physical be- 
havior. In the Kondo problem for instance, the crossover 
takes place between a weakly coupled two level system 
(the spin 1/2 impurity) at high energy, and a strongly 
coupled screened impurity at low energy, and is charac- 
terized by the Kondo temperature Tk |l| • This crossover 
has been widely studied in equilibrium It has also 
been considered in non equilibrium situations, both ex- 
perimentally [2j, |3[ and numerically [J, [5| . 

Many theoretical methods have been developed to at- 
tack quantum impurity problems analytically @, 0] , and 
it is reasonable to expect that progress can be made in 
the study of quench dynamics as well. Note however 
that the most potent method - the use of conformal in- 
variance [H - does not apply in the case of crossovers, 
since these precisely describe the departure from scale 
invariance. This means that in quantum impurity prob- 
lems, only the very long time, low energy, behavior can 
be described by methods such as those in || Q . 

Of all the quantities one may want to study in 
quenches, the most fundamental is probably the overlap 
of the states before and after the quench [10( . In equilib- 
rium, this is deeply related to the Anderson orthogonality 
catastrophe [Hj]. The type of situation we have in mind 



here is rather the work distribution in a situation where 
an impurity is suddenly coupled to an electron bath, the 
coupling being characterized by a crossover temperature 
7b in equilibrium. We expect that the work distribution 
should be a function of W/Tb, barring some scaling vi- 
olations such as those observed in the screening cloud 
problem [l2j |. 

We discuss here the main ideas of our approach, and 
illustrate them in the case of the resonant level model 
(RLM) . The Kondo model can be discussed similarly, at 
the price of some extra technical difficulties that will be 
presented elsewhere. We immediately warn the reader 
that, although the RLM is known as a 'non-interacting 
model', the quench is an extremely non-trivial operation 
in the free electron basis, as it affects an infinity of mul- 
tiple particle states. The technical difficulties involved in 
its study are entirely analogous to those occurring in the 
Kondo case. 

General quantum mechanics considerations show that 
the Loschmidt echo and the work distribution function 
can be expressed in terms of the scalar products of eigen- 
states of the Hamiltonian before and after the quench. 
The main point of our approach is that these scalar prod- 
ucts can all be determined analytically, in the limit where 
all energy scales are much smaller than the cutoff (band- 
width) , by using an axiomatic formulation very much like 
what is done in the form-factor approach to correlations 
in integrable massive quantum field theories. The cor- 
responding theoretical tools were put forward in earlier 
papers 
culties. 
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14j . There are however several major diffi- 
First, an infinity of such matrix elements has 
to be taken into account. Second, the associated sums 
over eigenstates are strongly divergent at low energies, 
via an infrared catastrophe inherent to the fact that we 
are dealing with a massless bulk theory. As we shall see, 
these difficulties can be controlled using renormalization 
tricks, and remarkably accurate results from the low to 
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FIG. 1: Quantum quench in the Resonant Level Model. 



the high energy regions can be obtained. 

Resonant Level Model. The spinless RLM involves 
two independent one-dimensional wires connected by 
tunneling through a quantum dot (the 'impurity'). The 
RLM is equivalent to the anisotropic Kondo system at 
the Toulouse point, and to the problem of an impurity in 
a Luttinger liquid with parameter g — ^, and therefore 
the approach outlined below will apply to these systems 
as well. After unfolding the wires, the Hamiltonian reads 
H (7) = J dxH(j), with the Hamiltonian density 



H(j) = -iv F i>\dxi>a 

a=l,2 



a) 

Here, the label a denotes the two wires, 7 is a tunneling 
amplitude (which we took, without loss of generality, to 
be the same for both wires). The equilibrium physics 
of the RLM is rather simple. It is convenient to define 
tf>± = "^("01 i ^2)1 so that "0- decouples from the impu- 
rity. The scattering matrix of the remaining fcrmion ip + 
on the impurity then reads S+{u)) — . The tunnel- 

ing term is a relevant interaction, thus creating an energy 
scale Tb = and the system flows under renormaliza- 
tion from the 7 = fixed point (independent wires) to 
a strong coupling fixed point 7 — 00 where the impurity 
is completely hybridized with the wires. At low energy, 
the only remaining effect of the impurity is a phase shift 
ip+(0 + ) = c 2il V + (Cr) with e 2iS = S+(lu = 0) so S = § . 

Quantum quench. The main goal of this paper is to 
understand analytically the quantum dynamics of this 
system after suddenly turning on (or off) the tunneling 
7. Let Ho = H(~f = 0) be the Hamiltonian of the system 
for t < 0, and Hi = H(j) the Hamiltonian for t > (see 
Fig. [TJ. The framework presented here is quite general 
and can be applied at finite temperature, but for sim- 
plicity, we will only consider the case T = and imagine 
that the system is initially prepared in the groundstate 
|«F(0)) = Vo 0) ) of H Q for t < 0. The wave function of 
the system at time t is then = e~ iHlt |*(0)). 

Work distribution and crossovers. We are interested 
in the work performed on the system during the quan- 
tum quench. In order to determine the work done W, two 
energy measurements are necessary, so that the work is 
not an observable but is rather characterized by a dis- 



(Ei — E U> )), where ^>o"i an d Eq 1 { are the eigenstates 
and energies of the Hamiltonian Hq before (resp. Hi af- 
ter) the quench. The work distribution P(W) has gath- 
ered a lot of attention recently [1, [^, [3, [l6| , and is 
especially of interest since it can be measured experi- 
mentally by spectroscopy (see also [13] )• The gen- 
erating function of the moments of the work G(t) — 
J dWe~ tWt P(W) is given by the Loschmidt echo 

G(t) - <^o°V Ho *e- lffl *l4 0) ): whi ch is, up to a phase, 
nothing but the overlap (^(^^(O)). For energies much 
smaller than the cutoff (bandwidth) , the Loschmidt echo 
and the work distribution take the universal forms G(t) = 
g(tT b ), and P{W) = ip(^), where SE = E^-E^ 
is the minimal work that should be performed on the 
system during the quench. They show some interesting 
crossover dynamics at t ~ T, (W — SE ~ Tb) (here, 
Tb is the equivalent of Tjc). This contrasts with the case 
of gapless quantum spin chains for which the Loschmidt 
echo at low energy is a pure power-law with an exponent 
related to the central charge [13, [20( • The time evolu- 
tion after the quench follows the Renormalization Group 
(RG) flow all the way from the UV (t < TJT 1 , weak 
coupling) to the IR low energy limit (t ^> T b ~ , strong 
coupling). The computation of a quantity such as G{t) 
thus requires all the non-perturbative physics in order 
to observe the crossover at t ~ Tr . Note that even in 
a free-fermion problem such as the RLM, the computa- 
tion of G(t) is highly non-trivial because it is non-local 
in terms of the fermions. 

Mapping onto a boundary problem. In order to com- 
pute the Loschmidt echo G(t), we first map the quantum 
impurity system onto an equivalent boundary problem. 
This can be done in general, and the Loschmidt echo 
G{t) can then be interpreted as the partition function of 
a 2D statistical problem, critical in the bulk, with non- 
conformally invariant boundary conditions. In our ex- 
ample of the RLM, after folding the -0 + -fermions system 
and decomposing ip + into Majorana (real) components, 
one finds that the RLM is equivalent to two indepen- 
dent copies of the critical Ising model, with a bound- 
ary magnetic field hb = This is consistent with 
the fact that the scattering matrix S+(ui) has the form 
of the reflection matrix in the Ising field theory [2lj |. 
The Loschmidt echo after a Wick rotation then reads 
G(t = —it) = (Zh b (r)) 2 , where Zh b (r) is the partition 
function of the critical Ising model in the half-plane, with 
finite boundary magnetic field hb for imaginary times be- 
tween and r, and free boundary conditions (BC) else- 
where (see Fig. [2]). Therefore, the Loschmidt echo follow- 
ing a quench in an integrable impurity problem reduces to 
the computation of a partition function in an integrable 
classical system with alternating zero and finite bound- 
ary magnetic fields. Note also that Z} lb (r) is normalized 
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FIG. 2: The Loschmidt echo in the resonant level model can 
be thought of as a partition in the Ising model with alter- 
nating vanishing and finite boundary magnetic fields Jit,. The 
quantum quench then effectively creates an arbitrary number 
of fermionic excitations in the system, r is the imaginary 
time. 



so that Z hb (r = 0) = Z hi =o(r) = 1. 

Low energy limit. The computation of Zh b (t) is non- 
trivial because the non-zero boundary magnetic field cre- 
ates an energy scale Tj = h%, and thus the problem 
cannot be solved using Boundary Conformal Field The- 
ory (BCFT) techniques. However, in the low-energy 
limit rh\ ^> 1, the BC becomes conformally invari- 
ant and Zh b (r) coincides with the two-point function 
of the Boundary Conditions Changing (BCC) operator 
from free to fixed BC in the Ising model. The scal- 
ing dimension of the latter is known to be /i = ^, 
so that Zh b (r) ~ (/i^t) -1 / 8 in that limit. Going back 
to the Loschmidt echo, this means that G(t) ~ i -1 / 4 
when tTf, 3> 1, where the exponent can also be inter- 
preted in terms of the phase shift of the ip + fermions 
2h = \ (^) = |. The vanishing of G(t) can therefore be 
traced back to the well-known Anderson orthogonality 
catastrophe. In terms of the work distribution, we thus 
obtain an edge singularity at low energy 0] 
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with a = (6/ it) — 1/4 in the RLM. This conclusion 
holds also for quenches in interacting impurity models, 
and one finds for example a = g/2 for the anisotropic 
Kondo case (5 = 1 being the isotropic model, and g =1/2 
the Toulouse point); and a — 1/8 for an impurity in a 
Luttingcr liquid, regardless of the Luttinger parameter. 

Form Factor approach. Although BCFT techniques 
can be used to analyze the low-energy limit of the quan- 
tum quench, the full non-perturbative computation of 



G(t) is much more involved. When hb is finite, one can 
still think of Zh b (r) as the two-point function of some 
generalization of BCC operators to non-conformal field 
theory. For integrable systems such as the ones we are 
considering, it should come as no surprise that the ma- 
trix elements of these operators can be handled [H, El 
using axiomatic techniques very similar to those used in 
the Form Factor (FF) approach to bulk correlations [13] • 
In particular, we find that the quench of the boundary 
magnetic fields creates an arbitrary number of fermions 
in the Ising field theory (see Fig. ^ , so that one obtains 
the formal expansion 



Zh b {r) 
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,/Wol 2 
(3) 

where ft is the rapidity of the i th fermionic excitation of 
, /3 n |fi) is the overlap between the 



energy e ft . hb (/3i 
groundstate of the Ising model |Q) with free boundary 
conditions, and . . . ,0 n ) h denote the eigenstates of 
the model in the presence of a finite boundary magnetic 
field hb- This inner product can be thought of as a FF 
of BCC operator in the boundary Ising model, which 
is known exactly from the axiomatic approach 3, 14 1. 
(For other works using integrability ideas in the context 
of quenches, see HHH]). 

Now, expression ^ involves an infinite sum of terms. 
Since the problem is massless in the bulk, n particle terms 
do not necessarily give rise to smaller contributions, as 
their energy can be arbitrarily low. The convergence of 
the expansion ([3]), or even its numerical usefulness, are 
thus not obvious a priori. Moreover, the integrals in- 
volved in this expansion are in fact IR divergent! It turns 
out nevertheless that the sum can be regularized by re- 
organizing the terms in the expansion. Details will be 
given elsewhere. We only quote the leading-order final 
result, 



\nZ hb (T) = 



du 1 
— 9x (uMu) + 



2! J 2tt Ui 



du 2 

2lTU2 



g x (u 1 +u 2 ) 



ui - u 2 

Ul + U 2 



- 1 *(« 1 )*(« a )+. . . . 



where x = rTb (recall r is the imaginary time) is an 
effective coupling parameter, dimensionless as expected, 
9x{u) = e~ xu — 1, and the Kernel ty(u) reads 



1 + u< 



■ exp 




cos^t 
cosh | sinh ~ 



(4) 

with ^(0) = 1 and "J (it) ~ u 1 as u — > oo. As a non- 
trivial check, this expansion can be resummed exactly 
in the IR limit x 3> 1, which yields Zh b (j) ~ x -1 / 8 
as expected from CFT. Moreover, for any value of x, 
Zh h (r) can be estimated by keeping only a few terms in 
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FIG. 3: (a) Loschmidt echo as a function of x = tTb, comparison between numerical results and Form Factors. Note that we 
rescaled the Loschmidt echo |G(£)| 2 in the main figure by a factor x 1 ' 2 in order to cancel the asymptotic power-law behavior 
for x = tTb ^ 1. (b) Universal work distribution from FF. At low energy, the work distribution has an edge singularity with 
an exponent a — (— ) 2 = i given by the dimension of a BCC operator while for W — SE ~ Tt, the system shows an interesting 
crossover physics. 



the FF expansion. Typically, the two and three-particle 
contributions are sufficient to obtain a very good approx- 
imation all over the crossover region. Finally, we note 
that, although we have used imaginary time to allow for 
a statistical mechanics interpretation of the results, ev- 
erything works as well in real time, with r = it. 

Lattice model. In order to validate the FF approach, 
we compare our analytical results to numerics in the 
RLM case |T]). We consider the itinerant fermion model 
(equivalent to a XX spin- 1/2 chain) with two weak links, 

H = -J^2(4 +1 a + h.c.)-(J'- J)(clc-i + c\cQ + h.c.), 

i 

(5) 

with J — 1 so that the Fermi velocity is vp — 2. At suffi- 
ciently low energies J' <C J = 1, the system is described 
by the effective Hamiltonian (JXJ) , with 7 cx J'. Given the 
free fermionic nature of the problem, it is even possible to 
identify exactly the energy scale T b ~ 2 J' 2 / J [26]. The 
Loschmidt echo G(t) following a quench from J' = to 
,/' 7^ can be expressed as a determinant that we eval- 
uate numerically - see [2(| for similar calculations. We 
compute G(t) on L = 1025 sites for different values of J', 
and find that the results indeed collapse onto a universal 
curve after rescaling the time scale by a factor 

The universal curve of the Loschmidt echo computed 
from numerics and from an extrapolated (see supplemen- 
tary material) two-particle FF expansion are shown in 
Fig- G3-(a). Note that because of finite size effects, one 
expects the curves for small values of J' to describe well 
the universal curve for small x only. We find that the FF 
expansion is in very good agreement with our numerical 
results, even in the interesting non-perturbative region 



t ~ J 1 " 1 where G(t) has a non-trivial behavior - note 
that there is no free parameter in the results, which must 
match without possible rescaling of the time axis. The 
resulting universal work distribution is shown in Fig. [3} 
(b). As expected, we observe an edge singularity © at 
low energy. Note that this power-law singularity would 
be smeared at finite temperature. For W ~ SE + the 
work distribution has a bump which we interpret as a 
signature of a Kondo resonance. 

Discussion. The new non-equilibrium dynamics ap- 
proach presented here is based on a generalization of the 
form-factor approach. The calculated time evolution of 
the Loschmidt echo was found to be in beautiful agree- 
ment with independent numerical studies of an equivalent 
lattice model. Applying the same formalism to the the 
Kondo case is a bit more involved, since the form factors 
are more complicated in this case. However, the problem 
is not fundamentally different from the case of the RLM, 
the physics of the quenches being essentially interacting 
in that case as well. 

Crossovers in quench dynamics are just beginning to 
get studied experimentally, especially in the context of 
quantum impurity problems. Recently, in Ref. Q a 
Kondo crossover scale in the absorption energy of a 
quenched quantum dot system was reported, reminiscent 
of our results in the time domain (Fig.[3]-(b)). One could 
well imagine analogous experiments in the context of cold 
atom systems, whereby a local quench is applied to the 
trapping potential, and the subsequent time evolution of 
the zero-momentum occupation number is monitored for 
associated crossover effects. 
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Supplementary material 

In this supplementary material, we provide more de- 
tails regarding the FF calculations. The problem is to 
compute the Ising matrix elements | h b (ft, ■ ■ ■ , fti|^) 1 2 
in the expansion 



n=0" 1=1 



x\ hb (A,...,^ n |n} | 2 . 



(6) 



Those matrix elements were computed exactly for an ar- 
bitrary change of BC from fixed boundary magnetic field 
h a to hb in IH 14 1 . In the case h a hb > 0, fermions are 
created by pairs, whereas the case h a hb < requires an 
additional particle to account for the presence of a do- 
main wall. The intermediate case we consider here is 
h a hb = (free BC corresponds to h a — 0), and it allows 
for an arbitrary number of fermions. We will start from 
the case h a hb > 0, the FF of the corresponding BCC 
operator read 



\G ba Wl,...^2n)\ 2 = 



I h b <fi|n>, | 2 



i-n si nh^-A|$ (ft )|2 ntanh : 



2 ft " ft 



(7) 



i<j 



where \£l) h is the vacuum of the theory with boundary 
magnetic field h, T a = /i 2 = e^ a , Tb = /i 2 = e' 9b , and $i 
is the rapidity of the i th particle. The function |$(ft)| 
has the following integral representation 



$(ftft,ft,) 



x exp 



1 



cosh 



+ OC 



2 4 
d tCQ3(^t)-COs(^t) 



4 cosh | sinh | 



(8) 



In order to describe free boundary conditions, we should 
take the limit f3 a — > — oo. We find 



(9) 



with it = e^, and ^(u) is the kernel used the main text 
ty(u) = — 5 exp 



1 




cos^i 



t cosh | sinh | 



(10) 



Using this function, the FF for our case h a hb — can be 
expressed as 



\GbWi 



,ft0| 2 



hb (/3 2 „,...,ft|») 
it 



n*sn 



(ii) 



where Mi = e". We also recall that |f2) is the vacuum 
of the theory for free BC. In terms of these FF, the ex- 
pansion ^ becomes 



Zh b (r) 



(12) 

If we plug (|lll) into this expression, we find that the 
integrals diverge at low energy m = e^ — > 0. However, 
the ratio can be regularized via a formal expansion of the 
inverse denominator and proper regrouping of the terms. 
This is similar to the trick used in [27]. We find the 
well-defined expression 



du 



dui 

2lTUi J 2ttu; 



2_ ( e ~ x(«l+-U2) 



1) 



Ml — M2 
Ml + M 2 



i r — r — r — 

3! J 27TM1 Jo 27TM 2 J 27TM 3 



(e 



-x(u 1 +U2+U3) 



Ml - M 2 



Ml + M2 7 \U\ + M 3 



Ml - M 3 



- 1 *(mi)*(m 2 ) 



U2 ~ U 3 
M2 + M 3 



Ml — M 2 
Ml + M 2 



\ I / 



U1 + U3J \u 2 + u 3 

I 



*(mi)*(m 2 )*(m 3 ) + ... (13) 



where we recall that x = rTb- Note that these IR thonality catastrophe as for h b — > 00, the inner product 
divergences can be traced back to the Anderson or- 
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oo(f2|f2)o should be exactly zero in the thermodynamic 
limit, so (jlll) becomes infinite in that limit. 

The regularized integrals in the FF are then computed 
numerically by keeping only a few terms (2 or 3 typically) 
in the expansion. The FF expansion is alternating, and 



averaging the expressions obtained at first and second 
orders has been shown to increase the precision by more 
than an order magnitude in the case where the analytical 
solution is known [27j . We obtain in this way an excellent 
approximation in the present case as well. 



